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ABSTRACT 


The  linear  theory  of  the  generation  of  electrostatic  bursts  of 
noise  by  electrons  trapped  in  chorus  wave  packets  is  developed  for  a 
finite  temperature  electron  beam  and  a  Maxwellian  electron  and  ion 
background.  The  growth  rates  determined  are  qualitatively  in  good 
agreement  with  those  obtained  by  previous  authors  from  a  more  idealized 
model.  Two  connected  instability  mechanisms  seem  to  be  occurring:  a 
beam  plasma  (electron-ion  two-stream)  instability  commonly  associated 
with  intensification  of  the  chorus  power  levels,  and  a  transitional  or 
borderline  resistive  medium  instability  commonly  associated  with  chorus 
hooks.  The  physical  reasons  for  the  two  mechanisms  is  discussed.  In 
the  second  case  electron  beams  are  difficult  to  identify  in  the 
particle  data.  An  expression  is  obtained  for  the  maximum  growth  rate 
in  terms  of  the  ratios  of  the  beam  and  electron  thermal  velocities  to 
the  beam  velocity,  and  of  the  beam  density  to  plasma  density.  It  is 
anticipated  that  this  may  allow  the  observed  peak  in  the  electrostatic 
noise  spectrum  to  be  used  as  a  diagnostic  for  the  beam  characteristics. 
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I.  INTRODUCTION 


In  recent  papers,  Reinleltner  et  al.  [1982;  1983]  presented  data 
showing  the  generation  of  electrostatic  waves  by  electrons  trapped  in 
chorus  wave  packets.  This  was  explained  on  the  basis  that  the  phase 
velocity  of  the  chorus  wave  packet  changes  with  the  local  plasma  con¬ 
ditions  and  magnetic  field,  resulting  in  acceleration  of  the  trapped 
particles.  The  acceleration  process  was  considered  in  more  detail  in  a 
more  recent  paper  [Gumett  and  Reinleitner,  1983],  Acceleration  shifts 
the  trapped  particles  in  velocity  space,  so  that  they  produce  a  "bump" 
(beam)  superimposed  on  the  untrapped  or  background  electrons.  Thu6 ,  a 
beam  instability  results  that  can  drive  electrostatic  waves. 

The  model  that  Reinleitner  et  al.  used  to  compute  the  spectrum 
produced  by  the  instability  was  that  of  a  resonance  (Lorentzian)  dis¬ 
tribution  of  electrons  and  a  delta  function  (zero  temperature)  beam 
distribution  in  velocity  space.  From  this  they  concluded  that  a 
"resistive  medium"  (negative  energy)  beam  instability  was  involved.  It 
is  the  purpose  of  the  present  paper  to  develop  the  theory  from  a  more 
realistic  model  with  the  objective  of  putting  the  interpretation  of  the 
data  on  a  firmer  foundation.  The  principal  conclusions  of  Reinleitner 
et  al.  [1983J  are  qualitatively  confirmed,  except  for  the  interpreta¬ 
tion  of  the  particular  instabilities  involved  in  the  electrostatic 


Vibi* 


bursts.  It  is  shown  that  the  electron-ion  two  stream  instability 
dominates  except  as  possibly  associated  wtih  the  chorus  hooks. 

The  theory  of  the  instabilities  is  developed  in  Section  II.  A 
discussion  of  the  predictions  and  the  data,  as  well  as  the  physical 
mechanisms,  is  given  in  Section  III.  The  results  are  summarized  in 


Section  IV 
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II.  THEORETICAL  MODEL 


One  example  of  the  electron  distribution  function  seen  when  chorus 
and  the  electrostatic  bursts  are  present  is  presented  in  Figure  1,  It 
consists  of  the  following  components: 

(1)  A  near  Maxwellian  central  body  with  a  thermal  velocity  ve  ~  5  x 
108  cm/sec.  (Note  that  it  does  not  fit  a  Lorentzian  well — see 

the  difference  in  the  log  plots  of  the  distribution  in  Figure  2.) 

(2)  Two  field  aligned  beams  with  beam  velocities  v^  ~  ±1.6  x  108 
cm/ sec. 

(3)  A  non-Maxwelllan  Call  which  sets  in  when  f(v)  is  less  than  lO-8 
times  its  value  at  the  peak,  and  at  a  velocity  over  twice  vj,. 

Since  the  non-Maxwellian  tail  makes  only  a  minor  contribution  to 
the  distribution  function,  it  should  be  quite  reasonable  to  neglect  it. 
Thus ,  we  may  take  the  reduced  distribution  function  of  the  background 
electrons  to  be 

f  (v)  -  )e^2/ve  (1) 

e  e 

where  v  is  the  velocity  along  the  field.  The  ions  will  also  be  taken 
to  be  Maxwellian.  The  beam  electrons  will  be  taken  to  be  a  square 


distribution: 
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fb(v) 


A 

0 


J  -  2e  ip/ me  <  v  <  J +  2e<j>/me 

otherwise 


(2) 


where  vj,  is  equal  to  the  phase  velocity  of  the  chorus  wave  w/k,  and  $ 
is  the  wave  potential.  If  we  take  the  width  of  the  beam  in  velocity 
space  to  be  2Av,  where  Av  is  the  beam  thermal  velocity,  then  A  « 
1/24*. 

The  electrostatic  beam  modes  that  can  be  present  are  primarily 
magnetic  field  aligned,  hence  they  will  be  taken  to  satisfy  the  one¬ 
dimensional  magnetic-field-free  electrostatic  dispersion  relation: 


where  <Upa2  -  cube2  =  nt^/ingCo  and  n^,  ne  and  nj,  denote 

the  background  ion  and  electron  and  electron  beam  densities,  respec¬ 
tively.  Here  Z(c)  is  the  familiar  plasma  dispersion  function.  The 
last  integral  is  straightforward  to  evaluate  from  Equation  2.  The  beam 
modes  in  Equation  3  satisfy  u^/k  ~  vj,  (where  w  ■  0^  +  iy),  and  for  the 
chorus  trapped  electron  beams  vj,2  »  vj2,  so  that  we  may  use  the  large 
argument  expansion  of  the  plasma  dispersion  function  for  the  ions 
(Fried  and  Conte,  1964]: 


i 
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C»1  _,2  i  3 

Z'(c)  ~  -2i/^?e  5  +—(1+  —  +...) 

?  2< 


Thus,  the  dispersion  relation  becomes 


(A) 


1 


2  2 


vi 


2 


(o)  -  kvb>2  ~  k2Av2 


(5) 


We  will  consider  the  solution  in  two  particular  cases  [Briggs,  1 964 J . 
(1)  Resistive  medium  instability  (ve  >  vj,  >  >4ne /m^  ve) 

Assume  ve2  »  vj,2  (i.e.,  ve  >  3  vj,)  to  simplify  the  problem.  Then 
the  small  argument  expansion  of  the  plasma  dispersion  function  may  be 
used 


C«1  _  2  o 

Z*U)  ~  -2A  i^e  4  -2(1  -  2c  +  ...) 


Then  Equation  5  becomes 


(6) 


,22  2 

2u)  o)  .  0)0) 

k  v  o)  k  v 

e  e 


2,,2  2 
-o)  /k  v. 


e  + 


0)0)  .  2,,  2  2 
(_E|)  e-«/kvl] 

k  vi 


2 
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(o)  -  kv^)2  -  k2Av2 


(7) 
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This  can  be  solved  in  a  straightforward  manner  to  give 


kv. 


,22  “be 
i  [k  4v  + 


g  +  ID 


1/2 


(8) 


where 


g 


2u> 


pe 


.  2  2 

k  v 

e 


2 

(!)  . 

_  _Ei_ 
.2  2 

k  vb 


(9) 


and 


oi  v 

e 


(10) 


The  neglect  of  the  ion  contribution  to  £>  is  justified  since  v^,2  »  v^2. 
Now  Equation  8  can  be  further  simplified  since  D2  <<  g2,  allowing  a 
separation  of  the  real  and  imaginary  parts: 


,  _  r.  2.  2  .  t>e 

“r  =  kvb +  lk  Av  +  “ — -t—2  2 — T~rrri 

(  1  +  2u>  /k  v  -  u>  ./k  v,  ) 

pe  e  pi  b 


±  /tt  ui?  (v^/v  3)(uj2  /uj2)e 

be  b  e  pe  r 


_  2,  2 

^  “vb/ve 


£  £  L  £  £  £ 

U)  V .  U)  .  3/2  ~  .  2  u>  V ,  U)  .  i/2 

1  ♦  2(-?)C-|)  -  (-f-)l  v  *  »r2!^ri(i  -  2  -f v  -  -§Mi 


The  top  signs  represent  the  unstable  mode  of  interest, 


orowth  rate 


yCoitf)  is  plotted  for  various  ratios  v^/Vg  and  a  couple  o.  .afferent 
ratios  n  =  u,be2/uipe2  *  n^/ng  in  Figures  3  and  4. 

It  is  clear  from  the  plots  that  the  frequency  at  which  the  growth 
rate  is  maximum  increases  with  increasing  v^/Vg  and  q.  This 
relationship  can  be  made  explicit  by  setting  3y/ Suij-  **  0  and  solving  for 
o*r.  If  we  denote  the  resulting  ojj-  by  ub,ax  (i.e.,  the  frequency  which 
maximizes  y),  the  result  is 


-(n  +  3=y- )  +  [(n  +  3*-)  +  8(4-)(n  +  2^)]' 


*i*r) 


Thus  iiimiv  is  a  function  of  three  ratios:  (Av/vj,),  (v^/ve)  and  q. 
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In  the  case  that  is  being  considered  here,  vb/ve  <  1,  the  beams 
are  very  difficult  to  separate  from  the  bulk  of  the  distribution  func¬ 
tion  since  they  occur  close  to  the  center  of  the  distribution  function. 
In  addition,  these  electron  beams  are  expected  to  be  highly  variable  in 
time  so  that  they  will  often  not  be  sampled  by  the  LEPEDEA  plasma 
instrument.  However,  Equation  13  combined  with  other  indirect  data  may 
be  adequate  to  determine  the  properties  of  such  electron  beams  in  the 
resistive  medium  case.  ve  can  be  estimated  from  the  LEPEDEA  data.  Av 
can  be  estimated  from  the  amplitude  of  the  potential  or  electric  field 
of  the  chorus  [see  Equation  2J.  u^,e  can  be  estimated  from  the 
continuum  radiation  cutoff  (which  is  just  another  measurement  of  the 
electron  density),  and  oyj  is  determined  by  the  maximum  in  the 
spectrum  of  the  electrostatic  bursts.  Since  Equation  13  is  not  very 
sensitive  to  n  it  can  be  used  to  predict  the  ueam  velocity  vb.  If,  on 
the  other  hand,  we  estimate  vb  from  the  chorus  phase  velocity,  which  is 
determined  from  the  whistler  dispersion  relation  and  the  local  electron 
plasma  density  and  magnetic  field,  the  spectral  maximum  may  provide  an 
estimate  of  n  and  hence  the  beam  density  nb.  Thus,  the  electrostatic 
spectrum  may  provide  a  diagnostic  for  one  of  the  beam  properties  in 
those  cases  for  which  the  beam  is  obscured  by  being  too  close  to  the 
center  in  the  LEPEDEA  data. 

In  Figure  5  we  have  plotted  (qnax/t0pe  as  a  function  of  vb/ve  for 
vb/ve  <  0*5  and  different  values  of  q.  The  plot  exhibits  the  monotone 
increase  of  with  vb/ve,  as  well  as  with  n,  although  the  latter  is 

minor  when  vb/ve  is  small.  The  plot  is  also  very  suggestive  that  ^ax 
approaches  oipe  as  vb/ve  approaches  1. 


i 
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2.  Electron-ion  two  stream  instability  (ve  <  v^,) 

This  is  the  case  exhibited  in  Figure  1,  and  is  just  the  well-known 
beam  plasma  Instability.  The  case  will  be  treated  here  as  a  review  for 
those  who  may  be  unfamiliar  with  this  beam  instability.  We  will  assume 
vj,2  >>  ve2  to  simplify  the  analysis.  Then  we  may  use  the  small 
argument  approximation,  Equation  4,  to  the  plasma  dispersion  function 
in  Equation  5  to  obtain  the  dispersion  relation: 


1  = 


2  ,22  2 
u)  „  k  v  uj  . 

-P£  n  +  1 _ £)  +  -Pi 

2  ^  2  2  ;  2 

0)  U)  w 


2 


2  2  2 
(u)  ~  kv  )  -  k  Av 

b 


-  2/rri  [  (• 


2  2  2 

pe)e-uj  /k  ve 


,3  3 

k  v 

e 


2 

u)u)  . 

r. _ Pi 

'•.3  3 

k  v. 
l 


2.,2  2 

-u>  /k  vj  ■ 


(14) 


This  can  be  rewritten  in  the  form 


(  oi  -  kv  )  2  =  k  2&v  2 

D 


2 

^e 


h(w)  +  ip 


(15) 


where 


U)  U)  U) 


(16) 


and  V  is  given  by  Equation  10. 
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A  close  Inspection  of  Equation  15  shows  that  for  Av2/vb2  <<  1  and 
0  «  nb/ne  <<  1  there  is  only  a  solution  if  u  -  kvb.  This  allows  one  to 
expand  h(u>): 


h(ui)  =  h(kvb)  +  (  jj  -  kvb)h'(kvb) 
Now  for  k2vb2  <<  u>p e 2  we  have 
|h(kvb) | 

jh'(kvb)]  ’  2w  »  lu)  "  kvb* 


(17) 


(18) 


Thus,  in  that  case  the  second  term  of  Equation  17  is  totally  negli¬ 
gible.  However,  for  kvb  ~  wpe ,  h(kvb)  <<  1  and  the  second  term  of 
Equation  17  is  not  negligible.  In  both  cases  h(w)  <  0  and  the  right- 
hand  side  of  Equation  15  is  negative  for  Av2/vb2  <  n,  so  that  an  insta¬ 
bility  exists  for  all  u>  to  wpe  provided  that  criterion  is  satisfied. 

For  kvb  >  iDpe ,  the  right-hand  side  is  positive,  so  only  a  stable  propa¬ 
gating  mode  with  normal  Landau  damping  occurs. 

For  k2vb2  <<  o>pe2  we  have  h(uj)  -  h(kvb)  and  |h(kvb)|  »  D,  so  that 

Wr  =  kvb  +  <9(  nZ7/ h 2 )  (19) 


y  =  k{nvb/[l  +  |(-|)]  -  Av2} 
Vb 


1/2 


(20) 
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For  kvb  ~  u)pe, 


h(o)>  =  h(u>pe)  +  (iu  -  (jjpe)h'(o)pe) 


*  r-  —  f. 


e  • 

2‘ 


(oj 


u)  ) 

pe 


pe 


(21) 


If  we  substitute  this  back  into  Equation  15  and  set  6  “  u >  -  (Dpe,  we 
obtain  the  cubic  equation: 


26“ 


3v' 


2v: 


-  i D)  ui  6 2  “  2k26v26  -  U)  (  n  u>2  +  ID  k2Av2)  =  0  (22) 

pe  pe  pe 


Equation  22  can  be  solved  by  perturbation  methods  to  give 


fl  i  Ve  f2)  1/3  Av 

T  peL  *  n  3v 2  2  2 

b  b 


(23) 


r/3,rwi/3 

"peL-^-) 


2  2 


'  1  - 


2x1/3  Av' 


3v/ 


(24) 


D  is  just  the  normal  (electron  Landau  damping,  given  by  Equation  10. 

Equation  23  gives  the  frequency  at  which  the  growth  rate  is  a 
maximum,  so  that  it  is  analogous  to  Equation  13.  This  is  plotted  as  a 
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function  of  v^/ve  for  various  n  in  Figure  6.  In  principal,  one  might 
use  the  theoretical  prediction  to  try  to  extract  information  on  the 
beam  from  the  spectrum.  But  in  this  case,  the  beam  characteristics 
show  up  quite  clearly  in  the  particle  data,  as  in  Figure  1.  Comparison 
of  Figure  6  with  Figure  5  shows  that  switches  from  an  increasing 

function  of  vj,/ve  to  a  decreasing  function  near  v^/ve  ~  1. 


.  .. 
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III.  DISCUSSION 


The  electron  distribution  function  shown  in  Figure  1  has  v^/ve  ~ 
3.  The  wave  spectrum  during  the  same  interval  is  shown  in  Figure  8. 
The  peak  in  the  intensity  occurs  near  f  =  6  Hz,  or  f/fpe  “  uj/ Wpe  - 
1.0.  This  spectrum  runs  from  about  0.9  fpe  up  to  1.1  fpe.  This  is  to 
be  compared  with  the  frequency  of  the  growth  rate  maximum  of  ~  0.97 
ujpe  predicted  by  Equation  23  for  the  beam  plasma  instability  if  we 
estimate  n  ~  10“  3  from  Figure  1.  Thus,  this  case  is  in  good  agreement 
with  predictions  of  the  beam-plasma  or  electron-ion  two-stream  insta¬ 
bility.  The  spectrum  of  electrostatic  noise  is  seen  to  be  associated 
primarily  with  an  intensification  of  the  chorus  waves  and  not  with  the 
chorus  hooks.  Another  example  of  this  case  of  bursts  associated  with 
chorus  band  intensification  was  shown  in  Figure  2  of  Reinleitner  et  al 
| 1983J .  The  spectrum  runs  up  to  fpe  with  a  peak  near  0.8  fpe,  consis¬ 
tent  with  the  beam  plasma  instability. 

Several  examples  of  electrostatic  bursts  associated  with  chorus 
hooks  were  given  in  Reinleitner  et  al.  (1983],  The  spectrum  in  this 
case  Is  different  than  for  the  chorus  Intensification  associated 
bursts.  The  spectrum  does  not  usually  extend  up  to  the  plasma  fre¬ 
quency,  and  often  peaks  near  0.5  fpe  or  less.  Furthermore,  the  bursts 
tend  to  be  short-lived,  normally  having  a  lifetime  of  a  second  or  less 
The  spectrum  is  consistent  with  a  resistive  medium  instability  or  an 
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instability  on  the  transition  between  resistive  medium  and  beam-plasma 
(see  Figures  3,  4,  and  5). 

The  reasons  for  the  two  different  types  of  bursts,  apparently 
associated  with  two  distinct  characteristic  chorus  waves,  can  be 
explained  as  follows.  When  the  chorus  hooks  are  present,  the  spectrum 
is  spreading  rapidly  to  higher  frequencies.  As  this  happens,  w/k  ~  Ju 
also  increases  rapidly,  trapping  and  accelerating  particles  in  the 
right  velocity  range.  In  this  case,  the  frequency  at  the  beginning  of 
the  hook  may  be  sufficiently  low  that  w/k  lies  more  on  the  main  part  of 
the  distribution  function  than  the  tail.  Since  the  phase  velocity 
changes  rapidly  the  wave  quickly  traps  and  accelerates  enough  particles 
to  create  a  significant  bump  in  velocity  space  before  moving  very  far 
in  velocity  space.  Thus,  v^  <  ve  when  the  instability  develops.  After 
uj  ceases  to  rise,  the  trapping  rate  drops  drastically,  the  particle 
"bump"  quickly  becomes  detrapped,  and  the  instability  ceases. 

In  the  case  of  the  electrostatic  noise  associated  with  intensifi¬ 
cations  of  the  chorus,  the  frequency  of  the  chorus  is  normally  2  or  3 
times  that  at  the  beginning  of  the  chorus  hooks.  Thus,  the  wave  phase 
velocities  may  lie  more  on  the  tail  of  the  distribution  function  at  the 
outset,  and  get  accelerated  more  slowly  by  gradients  in  the  Earth's 
magnetic  field.  The  number  of  trapped  particles  may  be  smaller  than 
that  of  the  case  of  the  chorus  hooks.  When  the  chorus  intensifies,  it 
can  pass  the  threshold  (growth  exceeds  damping)  of  the  beam  plasma 
instability.  While  the  trapped  beam  particles  become  untrapped  and 
relax  to  lower  energies  due  to  quasllinear  effects,  the  chorus  wave  is 


continually  trapping  and  accelerating  particles  to  regenerate  the  beam. 
Thus,  the  Instability  persists  as  long  as  the  chorus  does. 

The  results  of  this  paper  can  be  compared  to  those  of  Reinleltner 
et  al.  [1983]  In  which  a  model  of  a  Lorentzian  electron  distribution 
and  a  cold  beam  were  assumed.  First  of  all,  the  use  of  a  Lorentzian 
electron  distribution  causes  an  overestimation  of  the  growth  rate  for 
the  resistive  medium  instability.  Secondly,  the  neglect  of  beam  ther¬ 
mal  effects  overestimates  the  beam  plasma  instability  growth  rate, 
although  the  error  is  rather  small  near  the  peak  growth  rate.  Thirdly, 
the  position  of  the  peaks  compares  favorably  with  that  of  the  present 
model  for  the  resistive  medium  case,  but  not  as  well  for  the  beam- 
plasma  case.  The  reason  for  the  latter  discrepancy  is  that  a 
Lorentzian  distribution  does  not  properly  treat  the  thermal  effects  for 
frequencies  near  the  plasma  frequency.  Equation  14  should  be  the  more 
accurate  expression  of  the  peaks  for  the  resistive  medium  case  and 
Equation  23  a  useful  expression  for  the  peaks  of  the  beam  plasma  case. 
Overall,  however,  the  Reinleitner  et  al.  model  gives  a  good  qualitative 
estimate  of  the  principal  features  of  the  instability. 


IV.  SUMMARY  AND  CONCLUSIONS 


The  growth  rates  for  a  finite  temperature  beam  of  chorus  trapped 
electrons  in  a  Maxwellian  plasma  were  determined  for  both  the  resistive 
medium  (vj,  <  ve)  and  beam  plasma  (v^  >  ve)  instability  cases.  These 
results  are  qualitatively  in  good  agreement  with  those  of  Reinleitner  et 
al.  [1983],  except  in  the  conclusion  as  to  the  relative  importance  of  the 
two  mechanisms.  The  beam  plasma  (electron-ion  two  stream)  instability 
was  found  to  be  associated  with  the  long  time  electrostatic  bursts 
generated  during  the  gradual  intensifications  of  the  chorus.  A  border¬ 
line  resistive  medium  or  a  transitional  resistive-medium  beam-plasma 
instability  is  associated  with  the  chorus  hooks.  The  physical  reason  for 
the  apparently  different  mechanisms  in  the  two  cases  was  explained  in 
terms  of  the  Initial  phase  velocity  of  the  wave  and  the  acceleration 
mechanisms  in  each  case. 

For  the  case  of  the  chorus  hooks,  beams  are  normally  not  detectable 
in  the  particle  data.  The  maximum  growth  rate  was  determined  as  a  func¬ 
tion  of  the  ratios  of  beam  and  electron  thermal  velocity  to  beam  velocity 
and  beam  density  to  plasma  density  for  the  case  of  a  resistive  insta¬ 
bility.  This  may  allow  the  observed  peak  in  the  electrostatic  noise 
spectrum  to  be  used  as  a  diagnostic  for  beam  properties  in  this  case.  In 
particular,  everything  in  the  expression  for  the  maximum  growth  rate  is 
capable  of  being  estimated  except  possibly  the  beam  velocity  or  the  beam 
density.  Thus,  the  spectrum  may  be  used  as  a  diagnostic  for  beam  velo¬ 


city  or  density  in  such  cases 
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FIGURE  CAPTIONS 

Figure  1  Perspective  plot  of  the  electron  distribution  in  velocity 

space  on  August  10,  1979,  from  1353-1355  UT,  determined 
by  ISEE-1  LEPEDEA  data. 

Figure  2  Log  plots  of  the  Lorentzian  and  Maxwellian  distribution 

functions.  The  Maxwellian  curves  in  while  the  Lorentzian 
curves  out. 

Figure  3  Plot  of  the  growth  rate  y  as  a  function  of  v^,/ve  for  the 

negative  energy  or  resistive  medium  beam  instability,  for 
n  “  10“ 3 . 

Figure  4  Plot  of  the  growth  rate  y  as  a  function  of  vi,/ve  for  the 

negative  energy  or  resistive  medium  beam  instability,  for 
n  “  2  x  10“ 2 . 

Figure  5  Plot  of  the  frequency  of  maximum  growth  rate  i%aX  as  a 

function  of  v{,/ve  for  the  resistive  medium  instability. 

Figure  6  Plot  of  the  frequency  of  maximum  growth  rate  as  a 

function  of  vjj/ve  for  beam-plasma  instability. 
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Figure  7 


Wave  spectrum  from  the  SFR  instrument  during  the  time 
period  of  Figure  1. 


A-G83-527 


Figure  3 


-G83-526 


